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Abstract 

We report results of proper-time evolution analyses of B decays 
at Belle, based on a 29.1 fb _1 data sample recorded at the T(4S) 
resonance. These results include measurements of sin20i, the lifetime 
of i?-mesons, and the mass-difference between two mass eigenstates 
of the B°-B° system (Am^). Am^ is measured using three different 
methods: fully reconstructed hadronic modes, a fully reconstructed 
semi-leptonic mode, and a partially reconstructed Z?* ± 7r T mode. 
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1 Introduction 



The main goal of a 5-factory is to first observe and establish CP viola- 
tion in 5-decays and then to precisely measure the CP violation angles of 
the Unitarity Triangle, as well as the lengths of its sides. The aims is to 
either confirm the Kobayashi-Maskawa model 0] for the mechanism of CP 
violation or find evidence of physics beyond the Standard Model. Analysis 
of the proper-time evolution of B decays is essential for measuring mixing- 
induced CP violation, and is one of the unique features of an asymmetric 
e + e~ collider at the T(AS). By combining a large number of accumulated B 
events with a clean environment to study 5-decays, the asymmetric-collider 
experiments can now provide the most precise measurements that involve 
proper-decay time analyses, such as lifetimes and the oscillation frequency 
in B°-B° mixing, which is characterized by the mass-difference between the 
two mass eigenstates of the B°-B° system (AmJ. 

In this article, we report on the results of proper-time analyses of B decays 
at Belle: measurements of sin20 1 0j, the lifetime of 5-mesons, and Am^. 

2 KEKB Accelerator and Belle 

KEKB || is an asymmetric e + e~ collider of 3 km circumference, which con- 
sists of 8 GeV e~ and 3.5 GeV e + storage rings and an injection linear ac- 
celerator for them. It has one interaction point where the e + and e~ beams 
collide with a finite crossing angle of 22 mrad. The collider has been op- 
erated with peak beam currents of 1230 mA(e + ) and 780 mA(e~), giving a 
peak luminosity of 5.5xl0 33 /cm 2 /sec (as of the end of 2001, the target is 
10 34 /cm 2 /sec). Due to the energy asymmetry, the T(45) and its daughter 
Z?-pair are produced with /?7 =0.425 along the electron beam direction (z 
direction) in the laboratory frame. The average distance between the two 
decay vertices of B mesons is approximately 200 /im. A total integrated 
luminosity of 47 fb" 1 was accumulated during the period between October 
1999 and the end of 2001, of which 4.0 fb -1 was taken off-resonance and the 
rest was taken on the T(45). In the analyses reported here, we use a 29.1 
fb -1 data sample, which contains 31.3 million BB pairs, recorded at T(45) 
between January 2000 and July 2001. 

Belle is an international collaboration consisting of ~300 physicists from 
~50 institutes in 14 countries (U.S.A., Russia, Asia, Australia, and Europe). 
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The Belle detector is a general purpose large solid angle magnetic spec- 
trometer surrounding the interaction point. 

Charged particle tracking is done by a silicon vertex detector (SVD) and 
a central drift chamber (CDC). The SVD consists of three layers of double- 
sided silicon strip detectors (DSSD) at radii of 3.0, 4.5 and 6.0 cm. The CDC 
is a small-cell cylindrical drift chamber consisting of 50 layers of anode wires 
(18 stereo wire layers), covering 17° < 9[ ab < 150°. The CDC is operated with 
a i?e(50%)+C*2-ff6(50%) mixture. The CDC also provides measurements of 
the energy loss with a resolution of a(dE/dx)=6.9%. 

Particle identification is done by three detectors: dE/dx measurements 
in the CDC, time-of-flight counters (TOF) and aerogel Cherenkov counters 
(ACC). The TOF system consists of 128 plastic scintillators. The time resolu- 
tion is 95 psec (rms). The ACC consists of 1188 aerogel blocks with refractive 
indices of between 1.01 and 1.03, depending on the polar angle. By combin- 
ing information from these detectors, the efficiency for identification is 
about 90% and the ir fake rate is 6% with the requirement P(K/tt) > 0.6. 

The electromagnetic calorimeter (ECL) consists of 8736 CsI(Tl) crystal 
blocks, 16.1 radiation length thick, and covering the same angular region as 
CDC. Electron identification in Belle is based on a combination of dE/dx 
measurements in the CDC, the response of the ACC, and the position, shape 
and total energy (i.e. E/p) of its associated Csl shower. 

The outermost detector, for the measurement of /x ± and Kl (KLM), 
consists of 14 layers of iron (4.7 cm thick) absorbers alternating with resistive 
plate counters (RPC). 



3 Measurement of sin2</>i 

The Standard Model (SM) predicts a CP violating asymmetry in the time- 
dependent rates for the initial B° and B° decays to a common CP eigenstate, 
fcp- In the case where fcp — (cc)K°, the asymmetry is given by 



A(At) 



R(B° -> f C p) - R(B° -> f C p) 

R(B° - f C p) + R{B» - f C p) 

— £/ sin 20i sin(AmdAt), (1) 



where £/ is the CP-eigenvalue of fop, At is proper-time difference of two B 
decays, and <p\ is one of the three internal angles of the Unitarity Triangle, 
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denned as <j> x = vr - arg {^&) (= /3) [§■ 

\ cb c d J 

The measurement requires the reconstruction of B° — > f C p decays, the 
determination of the 6-flavor of the accompanying S-decay (/tag), a mea- 
surement of At, and a fit of the expected At distribution to the measured 
distribution using a maximum-likelihood method. 

We reconstruct -B decays to the following CP eigenstates ||: J/ipKs, 
tp{2S)K s , XciK s , VcK s for \ f = -1 and J/ipK L for = +1. The J/ip and 
i/)(2S) mesons are reconstructed via their decays to (£ = /i,e). The 

ij){2S) is also reconstructed via its J/ipir + ir~ decay, the xa via its 
decay, and the r\ c via its K + K~n° and Ks('n~ + 'n~~)K~'K + decays. We also use 
B° — > J /ipK*°{Ks^°) decays, where the final state is a mixture of even and 
odd CP. We find that the final state is primarily £/ = +1; the £j = — 1 
fraction is 0.19 ± 0.04(stat) ± 0.04(syst) 0. 

Except for B° — > J/ipK L , we identify 5 decays using the energy dif- 
ference A.E = i^g 118 — i?b™ s m and the beam-energy constrained mass Mb c = 

J (-Ebcam) 2 — (Pb 1s ) 2 > where * s ^ ne cms Deam energy, and E^ and p™ s 

are the cms energy and momentum of the B candidate. The B meson signal 
region is defined as 5.270 < M bc < 5.290 GeV/c 2 and a mode-dependent 
window around AE = 0. Table |l] lists the numbers of observed candidates 
(N eY ) and the background (iVbkgd) determined by extrapolating the rate in 
the non-signal AE vs. M^ c region into the signal region. 

Table 1: Numbers of observed events (N ev ) and the estimated background 
(iVbkgd) in the signal region for each fcp mode. 



Mode 


N cv 


-Wbkgd 


J/^(£+e-)K S (TT+TT-) 


457 


11.9 


J/ip(e + £-)K s (ir°Tr ) 


76 


9.4 


i/;(2S)(£ + £-)Ks(TT + n-) 


39 


1.2 


iP(2S)(J/ijn + TT-)Ks(TT + n-) 


46 


2.1 




24 


2.4 


7] c (K + K-TT )K S (Tr + TT-) 


23 


11.3 


Vc (K s K-7r+)K s (7T+7r-) 


41 


13.6 


J/^K*°(K S TT°) 


41 


6.7 


Sub-total 


747 


58.6 


J/^j(l+l-)K L 


569 


223 



Candidate B° — > J/ipKi decays are selected by requiring ECL and/or 
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KLM hit patterns that are consistent with the presence of a shower induced 
by a neutral hadron. The quantity p c ™ is calculated with the B° — > J / ipK^ 
two-body decay hypothesis. The distribution is fitted to the sum of 
signal and background distributions. There are 569 entries in the 0.2 < 
pf° s < 0.45(0.40) GeV/c signal region with KLM (ECL) clusters. 

Tracks that are not associated with a reconstructed CP eigenstate decay 
are used for flavor tagging. Initially, the 6-flavor determination is performed 
at the track level: high momentum leptons from b — > c£~V, lower momentum 
leptons from c — > sl + v, charged kaons and A baryons from b — > c — > s, high 
momentum pions that originate from decays of B° — > (nn) + , and slow 

pions from D* + — > D°ir + . We use the MC to determine a category-dependent 
variable that ranges from —1 for a reliably identified B° to +1 for a reliably 
identified B°, and depends on the tagging particle's charge, cms momentum, 
polar angle and particle-identification probability, as well as other kinematic 
and event shape quantities. The results from the separate track categories 
are then combined to take into account correlations in the case of multiple 
track-level tags. This stage determines two event-level parameters, q and r: 
q takes values q = +1(— 1) when the tag-side B meson is more likely to be a 
B° (B°); r is an event-by-event flavor-tagging dilution factor, which ranges 
from r = for no flavor discrimination to r = 1 for unambiguous flavor 
assignment. It is used only to sort data into six intervals of r, according to 
the flavor purity. 

The wrong-tag probabilities, wi (I = 1,6), are determined directly from 
the data for the six r intervals using fully reconstructed, self-tagged decays, as 
described in Sec. 5.1. The total effective tagging efficiency is J2i fi(^- — 2wi) 2 = 
0.270 ± 0.008(stat)^Q ^(syst), where fi is the fraction of the events in each 
r interval. 

The vertex positions for the fcp and / tag decays are reconstructed using 
tracks that have at least one three-dimensional coordinate determined from 
associated r-<f> and z hits in the same SVD layer along with one or more 
additional z hits in the other layers. Each vertex position is required to be 
consistent with the interaction point profile smeared in the r-0 plane by the 
B meson decay length. The fcp vertex is determined using lepton tracks 
from J/ip or i/)(2S) decays, or prompt tracks from i] c decays. The / tag vertex 
is determined from well reconstructed tracks not assigned to fcp- Tracks 
that form a K$ are not used. 

We determine sin2</>i by performing an unbinned maximum-likelihood 
fit of a probability density function (pdf) to the observed At distributions. 
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The pdf is a convolution of the theoretical distribution with a At resolution 
function, including the effect the of background and wrong-tag probabilities. 
The result is 

sin 20! = 0.99 ± 0.14(stat) ± 0.06(syst). 
The systematic error is dominated by the uncertainties in vertex reconstruc- 
tion and resolution. Figure [1] shows the observed At distributions for the 
q£f = +1 (solid points) and q£f = —1 (open points) event samples. In 




At (ps) 

Figure 1: At distributions for the events with g£j = +1 (solid points) and 
q£f = — 1 (open points). The results of the global fit (with sin20i = 0.99) 
are shown as solid and dashed curves, respectively. 

Fig. 0(a) we show the asymmetries for the combined data sample that are ob- 
tained by applying the fit to the events in each At bin separately. The smooth 
curve is the result of the global unbinned fit. Figures ^|(b) and (c) show 
the corresponding asymmetries for the (cc)Ks (£/ = ~ 1) and the J/ipKi 
(£/ = +1) modes separately. The fits give sin 20i values to be 0.84±0.17(stat) 
and 1.31 ± 0.23(stat), respectively, which are opposite, as expected. Fit- 
ting to the non-CP eigenstate self-tagged modes B° — > D ( -*^~n + , D*~p + , 
J/ipK*°(K + n-) and D*~i + u, yields 0.05 ± 0.04. As shown in Fig. g(d), no 
asymmetry is seen, as expected. 
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Figure 2: (a) The asymmetry obtained from separate fits to each At bin for 
the full data sample; the curve is the result of the global fit. The correspond- 
ing plots for the (b) (cc)K s (£/ = -1), (c) J/ipK L (£/ = +1), and (d) B° 
control samples are also shown. The curves are the results of the fit applied 
separately to the individual data samples. 
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4 Lifetime Measurements 



For lifetime measurements, B and B mesons are fully reconstructed in 
the following hadronic decay modes: B° — > D + tt~, D* + tt~, D* + p~, J/ipKs, 
J/^K*°(K-tt + ), and B~ ->■ _D 7r~, J/ipK'. The reconstruction of J/^X 
is made in the way quite similar to that used for the CP analysis described 
above. Neutral and charged D candidates are reconstructed in the follow- 
ing channels: D° — > K~n + , K~7T + 7r°, K~7t + 7t + tt~ , and D + —>■ K~7t + tt + . 
£>*+ _> ^0^+ can didates are formed by combining a D° candidate with a 
positively charged soft pion. To reduce the continuum background, a se- 
lection based on the ratio of the second to zeroth Fox- Wolfram moments H 
and the angle between the thrust axes of the reconstructed and associated 
B mesons is applied mode by mode. The vertices for fully reconstructed and 
associated B's are reconstructed using the same algorithm that is used in 
the CP analysis. We find 7863 B° and 12047 B~ events within the AE-M hc 
signal regions after the vertexing and selections are applied. 

We extract the B° and B~ lifetimes by performing an unbinned maximum- 
likelihood fit to the At distributions of B° and B~ simultaneously. The pdf 
is similar to the one used in the CP analysis. However, we have devel- 
oped a more elaborate resolution function for lifetime measurements. It is 
constructed as the convolution of four different contributions: the detector 
resolutions for the fully reconstructed B and associated B, the additional 
smearing from the associated B due to the inclusion of secondary tracks and 
the kinematic approximation that the B mesons are at rest in the cms. We 
allow for outliers with a large-width Gaussian to take into account any im- 
perfection of the resolution function. The fit includes 10 free parameters (7 
for resolution function and 3 for outliers) besides r B o and t b -. We obtain||: 

t B o = 1.554 ± 0.030(stat) ± 0.019(syst) ps, 

t b - = 1-695 ± 0.026(stat) ± 0.015(syst) ps, 

t b -/t b0 = 1.091 ± 0.023(stat) ± 0.014(syst). 
Figure ^ shows the distributions of At for B° and B~ events in the signal 
region with the fitted curves superimposed. The dominant sources of the 
systematic errors are the uncertainties in the background At shape and the 
modeling of resolution function. 
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Figure 3: At distributions of neutral (top) and charged (bottom) B meson 
pairs with fitted curves. The dashed lines represent the sum of the back- 
ground and the outlier component, and the dotted lines represent the outlier 
component. 



5 Am^ Measurements 

Am,j can be measured from the oscillation pattern of the At distributions of 
the mixed (same flavor, SF) and unmixed (opposite flavor, OF) B pairs: 

N SF oc e - |A * l/r B" [1 - cos(Am d At)] , (2) 
iVor oc e~ m/T B» [1 + cos(Am (i At)] . (3) 

We published a result using inclusive dileptons based on an early data 
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sample [ID]. Here, we report on the results using three other different sam- 
ples: fully reconstructed hadronic modes, a fully reconstructed semi-leptonic 
mode, and a partially reconstructed D* ± 7r T mode. 



5.1 Fully Reconstructed Modes 

We use fully reconstructed, self-tagged decay modes in two categories: hadronic 
decays (D l -*^n + , D*~ p + ) and semi-leptonic decay (B° — > D*~l + v). For 
hadronic modes, the same event samples are used as lifetime measurements. 
For B° — > D*~£ + v, we use the same D*~ decay chains as for the hadronic 
modes. The D*~ candidates are combined with /i + or e + candidates having 
charge opposite to that of the D*~ candidate. We exploit the massless char- 
acter of the neutrino using the missing-mass squared in the cms defined by 
ML = MM - C cos 6 b ,dh where MM = (E B - E DH f - \p B \ 2 - \pdh\ 2 , 
C = 2\pb\ \pdh\- Eb and ps [Ed*i and Pdh) are cms energy and momen- 
tum of B° (D* and lepton system). 9b, dh is an angle between vectors ps 
and Pd*i- It gives cos6 B) d*i = —MM/C by setting M^ iss = 0. We se- 
lect D*~£ + v candidates by requiring \ cos 6 B,D*t\ < 1-1- The 6-flavor of 
the accompanying B meson is assigned according to the flavor-tagging al- 
gorithm described above. The values of wi are obtained from the amplitudes 
of the time-dependent B°B° mixing oscillations: (iVoF — A^f) / (-^of+^Vsf) = 
(1 — 2wi) cos(ArrtrfAt). We perform a maximum-likelihood fit to SF and OF 
At distributions simultaneously with Arrid and w^s as free parameters. (Am^ 
is fixed at the world average value |nj when we determine Wi's.) We obtain 
Arrid = 0.521 ± 0.017(stat) iaoii(syst) ps -1 (prelim.) for hadronic modes 
and Arrid = 0.489 ± 0.012(stat) t^oi^syst) ps -1 (prelim.) for semi-leptonic 
mode. Dominant sources of systematic errors for hadronic modes are uncer- 
tainties in the resolution function and the background components, including 
mixing. Those for the semi-leptonic mode are uncertainties in the effects of 
the backgrounds. 

5.2 Partially Reconstructed D* ± tt =f 

In this method, B° — ► D^-n^ decays are reconstructed using only the fast 
pion (iTf) directly from the B° and the soft pion (tc s ) from D* + — ► D°tt + . 
Here, we use a method based on reconstructing the D° missing mass. We 
select a 7Tj with cms momentum 2.05 < p* < 2.45 GeV/c and an oppositely 
charged 7r s with a cms momentum of less than 0.45 GeV/c. The missing 
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mass is then calculated assuming the B meson is at rest in the cms. We 
select events with missing mass larger than 1.85 GeV/c 2 as D* ± 7r T candi- 
dates. We require R2 to be less than 0.5 to reduce the continuum back- 
ground. In order to identify the flavor of the other B, we require an addi- 
tional high momentum lepton with a cms momentum larger than 1.1 GeV/c. 
To reduce the background from incorrect tags (i.e. itf and leptons from the 
same B), we further require cos0| > —0.8 where 0| is the angle between 
a 717 and a lepton in the cms. We select 3686 OF (1213 SF) candidates 
with 73% (62%) purity after the missing mass cut and lepton-tag. The z- 
vertex of leptons is determined from the intersection of the lepton tracks 
with the profile of B\ decay vertices. An unbinned maximum likelihood fit 
is performed to the OF and SF Az distributions to extract Arrid- We ob- 
tain Arrid = 0.505 ± 0.017(stat) ± 0.020(syst) ps _1 (prelim.). The dominant 
sources of systematic errors are uncertainties in the background fraction and 
the resolution function. 

Figure ^ shows the OF-SF asymmetries as a function of Az for each 
method mentioned above together with the fit curves. 



6 Conclusion 

Belle (along with BaBar) presented the first significant observation of the 
CP violation in 2001 summer. The measured value is 

sin 20i = 0.99 ± 0.14(stst) ± 0.06(syst) 
with more than 6a significance. We have also measured the lifetimes of 
B° and B~ and Am^. The results are consistent with those of other ex- 
periments |12|(T^]. We find t b -/t^o is different from 1 with more than 



3a significance. By combining four different methods, we obtain Arrid = 
0.490 ± 0.010 ps _1 (prelim.). These are the most precise measurements so 
far. 
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Figure 4: Asymmetry as a function of the proper decay time difference for (a) 
hadronic modes, (b) semi-leptonic mode, (c) partially reconstructed B° — > 
D*n, and (d) dileptons. Fit curves are also shown. 
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